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ABSTRACT 


In  this  study,  the  Landau  and  cyclotron  growth  rates  of  whistler 
mode  waves  in  Earth's  bow  shock  are  calculated  using  electron  distribu¬ 
tion  functions  obtained  with  the  fast  plasma  experiment  on  ISEE  2. 

Three  electron  distribution  functions  measured  within  the  transition 
region  of  the  shock  are  analyzed.  An  important  feature  of  these  elec¬ 
tron  distribution  functions  is  the  presence  of  a  field-aligned  beam  with 

.  n  '<«  i 

a(Tj^  > anisotropy .  The  beam  velocity  vector  is  directed  toward  the 
magnetosheath.  The  calculations  show  that  the  electron  distribution 
functions  spontaneously  generate  whistler  mode  waves  with  plasma  rest 
frame  frequencies  between  about  .1  and  100  Hz.  The  wave  normal  angles 
of  the  generated  whistlers  range  from  0°  to  the  resonance  cone  angle. 
Electromagnetic  Landau  resonance  and/or  cyclotron  resonance  contribute 
to  wave  growth  over  the  range  of  velocity  distributions  observed.  The 
waves  that  are  generated  by  the  normal  cyclotron  resonance  have  wave 
vectors  directed  toward  the  solar  wind  while  those  that  are  generated  by 
the  Landau  and  the  anomalous  cyclotron  resonances  have  wave  vectors 
directed  toward  the  magnetosheath.  The  predictions  of  the  study  are  in 
qualitative  agreement  with  observations  of  whistler  mode  waves  near 


Earth's  bow  shock 


I.  INTRODUCTION 


Although  shock  waves  in  collisionless  plasmas  have  been  studied 
for  more  than  two  decades,  the  physics  of  collisionless  shocks  is  still 
not  well  understood.  Much  of  the  current  theoretical  effort  is  aimed 
at  formulating  models  that  self-consistently  describe  the  observed 
shock  structures  and  the  shock  dissipation  mechanisms.  In  a  collision¬ 
less  shock  wave  the  upstream  unshocked  plasma  is  transformed  into  the 
downstream  shocked  plasma  within  a  spatial  length  that  is  smaller  than 
the  mean  free  paths  of  the  plasma  particles.  The  region  separating  the 
upstream  and  downstream  plasmas  within  which  the  dissipation  occurs  is 
called  the  transition  region.  Observations  of  fluctuating  electric  and 
magnetic  fields  associated  with  shocks  in  the  laboratory  and  in  inter¬ 
planetary  space  indicate  that  plasma  wave  turbulence  plays  an  important 
role  in  electron  and  ion  heating  at  the  shock.  For  information  and 
additional  references  on  the  subject  of  both  natural  and  laboratory 
shock  waves,  the  interested  reader  can  refer  to  Tidman  and  Krall 
[1971 1,  Biskamp  [  1973] ,  Galeev  (1976),  Formisano  [1977],  and  Greenstadt 
and  Fredericks  [1979], 

Earth's  bow  shock,  the  shock  wave  formed  by  the  interaction  of  the 
solar  wind  with  Earth's  magnetic  field,  provides  an  excellent  labora¬ 
tory  for  the  study  of  collisionless  shock  waves.  The  characteristics 
of  the  bow  shock  are  highly  variable  because  the  parameters  that 
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characterize  the  upstream  solar  wind  vary  over  a  broad  range.  In  addi¬ 
tion,  because  spacecraft  instruments  have  dimensions  smaller  than 
characteristic  shock  scale  lengths,  the  shock  fine  structure  can  be 
studied  in  detail  if  high  time  resolution  instrumentation  is  used  or  if 
the  spacecraft  is  moving  slowly  with  respect  to  the  shock. 

Among  the  experiments  on  board  the  ISEE  1  and  2  satellites  are  a 
fast  plasma  experiment  [Bame  et  al.,  1978]  and  a  plasma  wave  experiment 
[Gurnett  et  al.,  1978].  The  fast  plasma  experiment  is  specifically 
designed  to  measure  distribution  functions  in  one  dimension  in  approxi¬ 
mately  100  milliseconds  and  distribution  functions  in  two  dimensions  in 
one  spin  period  (3  seconds).  In  particular,  the  time  resolution  of  the 
fast  plasma  experiment  is  sufficient  to  allow  the  measurement  of  dis¬ 
tinctive  features  of  electron  distribution  functions  within  the  shock 
transition  region  with  negligible  aliasing  of  time.  An  important  new 
feature  of  the  electron  distribution  functions  measured  by  the  fast 
plasma  experiment  within  the  shock  transition  region  is  a  field-aligned 
beam  moving  toward  the  magnetosheath  [Feldman  et  al.,  1982;  1983],  The 
beam  is  observed  for  both  quasi-parallel  and  quasi-perpendicular  tur¬ 
bulent  shocks.  Quasi-parallel  and  quasi-perpendicular  refer  to  the 
angle  between  the  shock  normal  and  the  ambient  magnetic  field,  0Bn.  For 
a  quasi-parallel  shock,  0gn  is  greater  than  0°  but  less  than  about  45° 
while  for  a  quasi-perpendicular  shock  0gn  is  greater  than  about  45°  but 
less  than  90°. 

The  plasma  wave  experiment  observes  many  types  of  plasma  turbu¬ 
lence  in  association  with  Earth's  bow  shock.  This  turbulence  is 
similar  to  the  plasma  waves  observed  near  the  bow  shock  by  previous 


instrumentation  [Fredericks  et  al.,  1968;  Rodriguez  and  Gurnett,  1975]. 
Among  the  different  types  of  electrostatic  and  electromagnetic  waves 
observed  is  whistler  mode  turbulence.  Whistler  mode  waves  are  electro¬ 
magnetic  waves  that  propagate  at  frequencies  below  either  the  electron 
gyrof requency ,  a)g,  or  the  electron  plasma  frequency,  o)p.  In  the  solar 
wind  and  bow  shock  regions,  u>g  is  always  much  less  than  ajp.  Typically, 
the  whistler  mode  noise  is  most  intense  within  the  shock  transition 
region  and  is  observed  to  extend  short  distances  into  the  upstream 
solar  wind  and  much  longer  distances  into  the  downstream  magnetosheath. 

The  purpose  of  this  study  is  to  demonstrate  that  the  electron 
beams  measured  by  the  ISEE  2  fast  plasma  experiment  in  the  transition 
region  of  the  bow  shock  can  generate  whistler  mode  instabilities.  The 
generation  mechanisms  are  identified  and  the  subsequent  propagation  of 
the  whistler  mode  waves  are  discussed. 

In  addition  to  data  obtained  with  the  fast  plasma  experiment  and 
the  plasma  wave  experiment,  the  study  also  makes  use  of  data  obtained 
with  the  ISEE  1  and  2  magnetic  field  experiment  [Russell,  1978], 
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II.  TRANSITION  REGION  ELECTRON  DISTRIBUTION  FUNCTIONS 

The  electron  distribution  functions  analyzed  in  this  study  were 
measured  by  the  fast  plasma  experiment  within  the  transition  region  of 
Earth's  bow  shock.  In  many  instances,  it  is  possible  to  fit  the  data 
obtained  by  the  fast  plasma  experiment  to  simple  analytic  functions. 
Because  the  measurement  of  and  fitting  procedure  for  these  electron  dis¬ 
tribution  functions  is  described  in  detail  in  the  recent  reports  of 
Feldman  et  al.  [1982,  1983],  only  a  brief  summary  will  be  given  here. 

In  Figure  1,  slices  of  electron  distribution  functions  measured 
along  the  direction  of  the  magnetic  field  on  December  13,  1977,  by  the 
I SEE  2  fast  plasma  experiment  are  shown.  This  figure  is  reproduced  from 
Feldman  et  al.  [1982],  The  figure  depicts  parallel  slices  in  the  solar 
wind  (SW),  the  shock  transition  region,  and  the  magnetosheath  (MS). 

Data  for  the  full  two-dimensional  electron  distribution  function 
obtained  with  the  fast  plasma  experiment  suggests  that  the  solar  wind 
electron  distribution  function  is  a  nearly  isotropic  Maxwellian  at  this 
time  and  that  the  magnetosheath  electron  distribution  function  is  a 
nearly  isotropic  modified  Lorentzian  at  this  time.  As  the  satellite 
enters  the  shock  transition  region,  field-aligned  electron  beams  are 
observed.  The  beams  are  moving  toward  the  magnetosheath  and  the  beam 
velocity  increases  with  increasing  shock  penetration.  This  figure 
supports  the  fact  that  the  electron  distribution  functions  measured  in 
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the  shock  transition  region  have  both  a  flat-topped  component  and  a  beam 
component . 

The  analytic  representation  of  the  electron  distribution  function 
in  the  plasma  rest  frame,  F(v^,  v ) ,  is  constructed  from  two  measure¬ 
ments  of  the  electron  distribution  function  in  planes  parallel  (v()  and 
perpendicular  (v^)  to  the  ambient  magnetic  field  separated  in  time  by 
about  0.75  sec.  Measurements  of  the  electron  distribution  function  in 
planes  other  than  those  parallel  and  perpendicular  to  the  magnetic  field 
confirm  that  the  fits  to  these  two  data  sets  yield,  for  the  purposes  of 
this  study,  a  sufficiently  accurate  representation  of  F(v^,  V||).  The 
simple  analytic  function  used  to  fit  the  data  is  a  sum  of  a  convected 
Maxwellian  and  a  modified  Lorentzian.  The  convected  Maxwellian  models 
the  electron  beam  and  the  modified  Lorentzian  models  the  flat-topped 
component  of  the  electron  distribution  function.  These  models  provide 
excellent  fits  to  both  the  parallel  and  the  perpendicular  slices  of  the 
distribution  function. 

Specifically,  F(v^,  V|()  is  given  by 

F(vi,v  || )  =  Fj  +  F2 


where 


F 


1 


(1  +  (■ 


vl)2  +  V1  R,  P/R 
7  '  J 


is  the  modified  Lorentzian  and 


2  2 

-m(v  -  v.)  -m  v 

F2  =  SeXp( - 2^T - )exP(2^TT) 


is  the  convected  Maxwellian.  In  these  equations,  P  and  R  are  real 

pi ' 

numbers,  and  TB  are  perpendicular  and  parallel  temperatures,  m  is 
the  electron  mass,  k  is  Boltzman’s  constant,  and  Cj,  C2,  vj,  and  V3 
are  constants. 

It  should  be  noted  that  for  the  electron  distribution  functions 
analyzed  in  this  study,  the  direction  of  positive  vj,  be  it  parallel  or 
anti-parallel  to  the  magnetic  field,  is  determined  by  the  bow  shock 
geometry  at  the  observation  time.  The  positive  vN  direction  is  defined 
to  be  such  that  vfi»n  is  positive,  where  v H  is  a  unit  vector  in  the 
positive  V|  direction  and  n  is  a  unit  vector  in  the  direction  of  the 
shock  normal  and  toward  the  solar  wind.  Consequently,  because  the  beam 
velocity  vectors  are  always  directed  toward  the  magnetosheath,  the 
convection  speed,  V3,  is  always  negative. 

The  magnitude  of  the  beam  convection  speed  increases  with  increas¬ 
ing  shock  penetration.  This  variation  is  consistent  with  the  inter¬ 
pretation  that  the  beams  represent  solar  wind  electrons  that  are 
accelerated  by  a  macroscopic  electric  field  within  the  shock  transition 
region.  The  parallel  temperature  of  the  beam,  T g ,  is  obtained  from  the 
fit  to  the  fast  plasma  experiment  data.  However,  in  many  cases  the 
corresponding  perpendicular  beam  temperature,  T|,  changes  significantly 
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during  the  time  it  takes  for  the  detector  to  rotate  from  the  parallel 
to  the  perpendicular  orientations.  Consequently,  is  estimated 
assuming  that  the  first  adiabatic  invariant  of  the  electrons  is  con¬ 
served  as  the  electrons  move  from  the  solar  wind  into  the  shock.  When 
the  relative  motion  between  the  spacecraft  and  the  bow  shock  suffi¬ 
ciently  slow  so  that  Tj_  can  be  obtained,  the  calculated  and  n  ured 

perpendicular  beam  temperatures  are  in  good  agreement . 

In  this  study,  three  electron  distribution  functions  mea  a  by 
the  ISEE  2  fast  plasma  experiment  within  the  transition  region  of 
Earth's  bow  shock  are  analyzed.  Of  the  three,  two  were  measured  on 
December  13,  1977,  and  one  was  measured  on  October  15,  1978.  The  two 
electron  distribution  functions  obtained  in  the  transition,  region 
on  December  13,  1977,  correspond  closely  to  the  distributions  measured 
at  17:34:55.8  and  17:34:58.6  that  are  depicted  in  Figure  1.  The 
October  15,  1978,  electron  distribution  function  is  included  in  this 
study  because  the  magnitude  of  the  convection  speed  of  the  Maxwellian 
modelling  the  electron  beam  is  particularly  small,  which  leads  to  some¬ 
what  different  instability  conditions  than  for  the  December  13  cases. 

In  Figure  2,  contour  plots  of  the  best-fitting  analytic  models  for 
the  three  electron  distribution  functions  analyzed  in  this  study  are 
shown.  Plotted  are  contours  of  constant  logjg  F(vL,  vll)  spanning  one 
order  of  magnitude  of  F.  The  convected  Maxwellians  with  T ^  >  T n 
anisotropies  that  model  the  electron  beams  are  clearly  visible  as  is 
the  flat-topped  feature  of  the  modified  Lorentzians.  In  Table  1,  the 
constants  used  to  specify  these  three  electron  distribution  functions 
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are  given.  In  Table  2,  the  bow  shock  parameters  at  the  three  observa¬ 
tion  times  are  summarized.  For  both  shock  crossings,  the  bow  shock 
classification  is  turbulent. 


III.  METHOD  OF  CALCULATING  THE  GROWTH  RATE 


To  calculate  the  growth  rate  of  whistler  mode  waves,  the  results 
obtained  by  Kennel  and  Wong  [1967]  using  the  linearized  Vlasov  and 
Maxwell  equations  are  employed.  The  expression  used  is  equivalent  to 
the  results  given  by  Kennel  [1966].  If  Landau  and  first  order  cyclo¬ 
tron  resonance  terms  are  included  and  ion  motion  is  neglected,  the 
growth  rate,  is  given  by 


= 


'i  1  r 


“-|slg"<k7>  Hr  Te  I  vI(1-i  +  L  +  Vduj 

II  oi  ° 

r 


(1) 


where 


I  = 
m 


+  £rVi(x)  -  /7<77,£lJm(x))2(v!  3F 


1  9v , 


WU>g  3F_) 

Si  9V1 


W 


a)  +  mw 

=  _£ _ 1 

k 


and  W,  the  total  wave  energy,  is  given  by 


W(£,u>)  =  [£*(£,  a>)4(£,a,)  +  a>Kh)-£(k+,w)] 
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In  these  equations,  F  is  the  electron  distribution  function  normalized 

to  unity,  o)r  is  the  real  frequency  of  the  wave,  k  is  the  wave  number, 

e  ,  eD,  and  e„  are  the  left-hana,  right-hand,  and  parallel  electric 

field  polarization  vectors  of  the  wave,  K^1  is  the  Hermitian  part  of  the 

dielectric  tensor,  is  the  absolute  value  of  the  electron  gyro- 

f^i 

frequency,  and  x  is  - .  The  Landau  contribution  to  the  growth  rate 

U)g 

is  the  m  =  0  term  and  the  normal  and  anomalous  cyclotron  resonance  con¬ 
tributions  are  given  by  the  m  =  -1  and  the  m  =  +1  terms.  A  positive 
growth  rate  indicates  exponential  growth  and  a  negative  growth  rate 
indicates  exponential  damping. 

In  order  to  employ  Equation  1  to  calculate  growth  rates  for  whis¬ 
tler  mode  waves,  the  real  part  of  the  whistler  mode  dispersion  relation 
must  be  solved  for  ajj.  as  a  function  of  k.  For  simplicity,  this  study 
assumes  that  the  real  part  of  the  whistler  mode  dispersion  relation  is 
that  of  a  cold  plasma.  When  the  wave  phase  velocity  is  much  greater 
than  the  electron  thermal  speed,  hot  plasma  effects  on  dispersion  for 
the  whistler  mode  are  small  (e.g.,  see  Leimohn  [1962]).  However,  in 
the  work  reported  here  the  wave  phase  velocities  are  often  comparable 
to  the  electron  thermal  velocities.  To  investigate  the  error  intro¬ 
duced  by  this  simplifying  assumption,  the  numerical  code  of  Forslund 
[1979],  which  solves  the  fully  electromagnetic  linear  dispersion  rela¬ 
tion,  has  been  utilized  to  determine  the  real  frequency,  wr,  for  a 
given  wave  number,  k  [S.  P.  Gary,  private  communication].  For  para¬ 
meters  typical  of  this  study,  the  results  indicate  that  for  a  given 
wave  number  and  for  wave  normal  angles  0  <  45°,  where  the  wave  normal 


angle  is  the  angle  between  the  wave  vector  and  the  magnetic  field,  the 
real  frequency  for  a  cold  plasma  and  the  real  frequency  obtained  using 
the  code  of  Forslund  differ  by  not  more  than  30  percent.  If  k  is  held 


constant  and  ojj-  is  allowed  to  vary  by  30  percent,  the  computations  show 
that  the  change  in  the  growth  rate  calculated  using  Equation  1  is,  for 
the  majority  of  unstable  wave  numbers,  small.  As  the  wave  normal  angle 

0Jr 

approaches  the  resonance  cone  angle,  9r6S  =>  cos  1  ( — ),  the  index  of 

0)g 

refraction  becomes  infinite  and  the  wave  becomes  predominantly  electro¬ 
static.  For  wave  normal  angles  near  the  resonance  cone  angle,  the 
difference  between  the  real  frequency  for  a  cold  plasma  and  the  real 
frequency  obtained  using  the  code  of  Forslund  for  a  given  k  is  larger 
than  it  is  for  small  wave  normal  angles.  Consequently,  it  should  be 
noted  that  while  the  primary  results  of  this  study  are  not  affected  by 
the  simplifying  assumption  used  in  calculating  the  growth  rate,  the 
calculated  growth  rates  for  wave  normal  angles  near  0Res  are  not 


accurate 


IV.  GROWTH  RATE  CALCULATIONS  AND  GENERATION  MECHANISMS 


Growth  rates  for  whistler  mode  waves  have  been  calculated  using 
Equation  1  and  the  electron  distribution  functions  depicted  in  Figure 
2.  The  results  indicate  that  these  electron  distribution  functions  are 
unstable  to  whistler  mode  waves.  In  this  section,  the  growth 
mechanisms  that  generate  the  whistler  mode  waves  are  discussed. 

From  Equation  1  it  is  observed  that  the  total  growth  rate  is  a  sum 
of  three  terms:  the  two  cyclotron  resonance  contributions  and  the 
Landau  resonance  contribution.  The  sign  of  the  growth  rate  is  positive 
if  the  sum  of  the  three  contributions  is  positive.  The  signs  of  the 
m  *  -1,  0,  and  +1  contributions  are  positive  or  negative  depending  on 
whether  or  not  the  following  quantity  is  positive  or  negative: 


2  , 

VldVi 


V»  = 


w  +  mo) 

r  _ 8 


(2) 


which  can  be  rewritten  in  the  form 
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Integrating  the  second  term  in  Equation  3  by  parts  and  introducing  the 
electron  pitch  angle,  a,  where  a  is  measured  from  the  positive  vt|  axis. 
Equation  3  can  be  rewritten  as 


u)  2to  oo 

sign  Or—)  f-  j  Fv  dv 
*»  kii  °  1  1 


T  3F  2.  , 

J  r-v.dv 
A  da  1  i 


to  +  mu) 


(A) 


The  whistler  mode  waves  generated  in  the  transition  region  by  the 
electron  distribution  functions  shown  in  Figure  2  fall  naturally  into 
two  categories:  waves  with  wave  vectors  predominant ly  in  the  opposite 
direction  of  the  beam  velocity  and  waves  with  wave  vectors  predomi¬ 
nantly  in  the  direction  of  the  beam  velocity.  For  waves  propagating  in 
the  direction  opposite  to  the  beam,  i.e.,  toward  the  solar  wind,  k(,  is 
positive  and  the  growth  mechanism  is  the  normal  m  =  -1  cyclotron  reso¬ 
nance.  All  three  electron  distribution  functions  shown  in  Figure  2 
generate  whistler  mode  waves  via  the  m  =  -1  cyclotron  resonance.  For 
waves  propagating  in  the  direction  of  the  beam,  i.e.,  toward  the  mag¬ 
netosheath,  k|j  is  negative  and  the  growth  mechanisms  are  the  m  =  0 
Landau  and  the  anomalous  m  =  +1  cyclotron  resonance.  Of  the  three 
electron  distribution  functions  shown  in  Figure  2,  only  the 
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distribution  function  measured  on  October  15,  1978,  generates  whistler 
mode  waves  propagating  in  the  direction  of  the  beam.  In  the  following 
subsections,  the  cases  k|  >  0  and  kN  <  0  are  examined. 

A,  Wave  Vectors  Directed  Toward  the  Solar  Wind: 

Generation  Via  the  m  =  -1  Resonance 

To  investigate  the  generation  of  whistler  mode  waves  by  the  normal 
m  *  -1  cyclotron  resonance  contribution  to  the  growth  rate,  consider 
Equations  2  through  4  when  m  =  -1.  Because  the  beams  depicted  in  Figure 
2  are  moving  in  the  direction  of  the  magnetosheath,  the  parallel  compo¬ 
nent  of  the  wave  vector  must  be  positive  for  interactions  via  the  m  =  -1 
cyclotron  resonance.  For  the  distribution  functions  shown  in  Figure  2, 

the  second  term  in  Equation  2  is  negative.  Therefore,  in  order  for 

3  3F 

Equation  2  to  be  positive,  v,  — —  evaluated  at  the  m  =  -1  cyclotron 

1  3v  || 

0)  -  (i) 

resonance  velocity,  vfl  *  — ^  ,  and  integrated  over  all  v^  must  be 

II 

sufficiently  positive  to  overcome  the  negative  term.  Equivalently,  Equa¬ 
tion  4  expresses  the  result  that  a  region  of  sufficiently  large 
3F 

anisotropy,  —  <  0,  along  the  m  =  -1  cyclotron  resonance  velocity  must 
3a 

exist  if  the  growth  rate  is  to  be  positive. 

From  these  facts  it  is  expected  that  the  electron  distribution 

functions  shown  in  Figure  2  will  generate  whistler  mode  waves  if  the 

cyclotron  resonance  velocity  has  a  larger  magnitude  than  the  magnitude 

of  the  convection  speed  of  the  Maxwellian  modelling  the  electron  beam, 

| v 3 1 .  In  this  region  of  the  electron  distribution  functions  shown  in 
9F  3  K 

Figure  2,  -r—  is  positive  and  —  is  negative. 


The  left  column  of  Figure  3  shows  representative  growth  rates 
obtained  using  the  electron  distribution  function  measured  within  the 
shock  transition  region  at  1734:55.7  on  December  13,  1977.  At  this 

0)e  Wjj 

time,  ■=*■  is  about  315  Hz  and  -s*-  is  about  20  kHz.  Plotted  are  the 
ratios  of  total  growth  rate,  u>i,  to  real  frequency,  u^,  as  a  function  of 
the  wave  normal  angle.  The  growth  rates  for  the  two  frequencies 
0JJ./21T  ■  110,  and  120  Hz  are  given.  This  figure  depicts  whistler  mode 
waves  with  positive  growth  rates.  The  absence  of  a  plotted  growth  rate 
indicates  damping.  For  the  frequencies  shown  here,  both  the  m  =  +1 
cyclotron  resonance  contribution  to  u>£  and  the  m  =  0  Landau  resonance 
contribution  to  are  negative  and  small  compared  to  the  positive 
contribution  of  the  m  =  -1  cyclotron  resonance.  The  left  column  of 
Figure  3  indicates  that  as  the  frequency  is  increased  the  range  of 
unstable  wave  normal  angles  decreases  while  the  centroids  of  these 
angles,  which  range  from  0  to  about  80  degrees,  increase.  This  trend 
can  be  understood  from  a  plot  of  the  m  =  -1  cyclotron  resonance 
velocity. 

The  right  column  of  Figure  3  shows  the  m  =  -1  cyclotron  resonance 
velocity  as  a  function  of  wave  normal  angle  for  the  two  frequencies 

a»r 

■s—  =  110  and  120  Hz.  For  the  electron  distribution  function  measured 
2tt 

at  1734:55.7  on  December  13,  1977,  the  Maxwellian  modelling  the  electron 
beam  is  convected  into  the  shock  at  a  speed  of  about  3255  km/s.  It  is 
expected  that  the  resonance  velocity  must  have  a  larger  negative  value 
than  -3255  km/s  if  instability  is  to  occur  as  only  then  will  9F/9a  be 
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sufficiently  negative  (or  3F / 3v  H  sufficiently  positive)  along  the  inte¬ 
gration  path.  This  point  is  illustrated  in  Figure  4  where  the  reso- 

<i>j- 

nance  velocity  range  giving  rise  to  instability  for  the  case  =  110  Hz 

C  TT 

as  obtained  from  Figure  3,  is  overlayed  on  a  contour  plot  of  the  elec¬ 
tron  distribution  function  measured  at  1734:53.7  on  December  13,  1977. 

3F 

The  resonance  velocity  lies  in  a  region  of  positive  — —  and  large 

dV 

3F 

anisotropy,  —  <  0. 

dot 


B.  Wave  Vectors  Directed  Toward  the  Magnetosheath: 
Generation  Via  the  m  =  0  and  +1  Resonances 


To  investigate  the  generation  of  whistler  mode  waves  by  the  m  =  0 
Landau  and  the  m  =  +1  anomalous  cyclotron  resonances,  Equations  2 
through  4  are  again  examined.  Because  the  convection  speeds  for  the 
electron  beams  shown  in  Figure  2  are  negative,  for  the  generation  of 


whistler  mode  waves  via  the  m  =  0  and  the  m  *  +1  interactions  kg  is  nega 

tive,  corresponding  to  wave  vectors  directed  toward  the  magnetosheath. 

From  Equation  2  it  is  seen  that  in  order  for  the  m  =  0  Landau  con- 

3  3F 

tribution  to  to  be  positive,  -v  — —  evaluated  at  the  Landau  reso- 

^  v  || 

U>j-  " 

nance  velocity,  vB  =  -r — ,  and  integrated  over  all  Vi  must  be  positive. 

II  k 

Taking  into  account  the  sign  of  kj  with  the  expression  — p,  this  will 

k  k  •jf  11 

occur  when  regions  of  large  positive  p^  j  exist  along  the 

integration  path.  Considering  the  m  =  +1  cyclotron  term  for  the  case 


when  kj  is  negative,  the  first  term  in  Equation  4  is  negative.  Conse¬ 
quently,  the  m  =  +1  cyclotron  resonance  contribution  to  the  growth  rate 
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will  be  be  positive  when  regions  of  large  anisotropy,  —  >  0,  exist 

OOt 

along  the  integration  path. 

Of  the  three  electron  distribution  functions  shown  in  Figure  2, 

only  the  distribution  function  measured  on  October  15,  1978,  generates 

whistler  mode  waves  via  the  m  =  0  and  the  m  =  +1  resonances.  This  is 

because  the  convection  speeds  for  the  electron  beams  in  the  distribution 

functions  measured  on  December  13,  1977,  are  large.  Because  the  convec- 

3F 

tion  speeds  are  large,  — —  is  near  zero  at  typical  Landau  resonance 

3vn 

velocities  and  the  Landau  contribution  to  the  growth  rate  is  small. 
Consequently,  the  m  *  -1  cyclotron  contribution  to  damping  dominates  any 
positive  contribution  to  the  growth  rate  due  to  the  m  =  +1  cyclotron 
resonance.  This  is  not  the  case  for  the  electron  distribution  function 
measured  on  October  15,  1978,  which  has  a  smaller  beam  velocity. 

Because  the  beam  velocity  is  small,  the  m  *  0  Landau  contribution  to 
is  large  and  positive. 

The  left  column  of  Figure  5  presents  a  sample  of  the  growth  rates 
obtained  using  the  October  15,  1978,  electron  distribution  function. 

u>r  Wg 

The  results  are  given  for  one  frequency,  •=—  *  80  Hz.  At  this  time, 

L  7T  L  TT 

is  about  110  Hz  and  is  about  20  kHz.  The  top  left  panel  shows  the 

o>i 

m  ■  0  Landau  resonance  contribution  to  the  ratio  —  while  the  bottom 

a>r 

U>i 

left  panel  gives  the  m  *=  +1  cyclotron  resonance  contribution  to  — . 

wr 

u>i 

The  m  ■  -1  cyclotron  resonance  contribution  to  —  is  negative  and 

negligible  for  these  wave  normal  angles  and  this  frequency.  The  right 

column  of  Figure  5  depicts  the  Landau  resonance  velocity  and  the  m  =  +1 

cyclotron  resonance  velocity  as  a  function  of  wave  normal  angle  for 
o>r 

tt-  -  80  Hz. 

2ir 
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From  Figure  5  it  is  observed  that  the  Landau  contribution  to  is 
a  maximum  at  an  intermediate  wave  normal  angle  and  decreases  to  zero  at 
0°  and  0Res  (0Res  “  43°  here).  The  contribution  is  positive  over  a 
wide  range  of  wave  normal  angles  because  the  Landau  resonance  velocity 

k„  3F 

lies  in  a  region  of  positive  -n — r  —  over  a  wide  range  of  wave  normal 

TV  3vh 

angles.  The  Landau  growth  rate  goes  to  zero  at  0°  because  the  whistler 

mode  wave  is  purely  transverse  at  that  angle  and  goes  to  zero  at  0Res 

because  the  index  of  refraction  for  the  whistler  mode  becomes  very 

large  near  0Res,  thereby  moving  the  resonance  velocity  out  of  the 

k,  3F 

region  of  large  positive  -rr— r  ~ — . 

'  k  II '  3v  II 

The  m  =  +1  cyclotron  resonance  contribution  to  is  positive  for 

wave  normal  angles  between  about  10°  and  42°.  The  m  =  +1  contribution 

falls  to  zero  as  the  wave  normal  angle  approaches  0°  because  I+j  in 

Equation  1  is  identically  zero  at  0°.  This  corresponds  to  the  fact  that 

a  parallel  propagating  whistler  mode  wave  is  right-hand  polarized  with 

respect  to  the  magnetic  field,  whereas  the  anomalous  m  =  +1  resonance 

requires  a  left-hand  polarization  component  for  a  resonant  interaction  to 

occur.  The  m  *  +1  contribution  is  zero  for  wave  normal  angles  greater 

than  about  42°  because  the  m  =  +1  cyclotron  resonance  velocity  does  not 

9F 

lie  in  a  region  of  sufficiently  positive  for  wave  normal  angles 
greater  than  42°. 

In  Figure  6,  the  range  of  unstable  resonance  velocities  for  the 
m  =  0  and  the  m  =  +1  resonances,  as  obtained  from  Figure  5,  is  overlayed 
on  a  contour  plot  of  the  October  15,  1978,  electron  distribution  func- 


t  ion . 


The  Landau  resonance  velocity  lies  in  a  region  of  positive 
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— —  while  the  m  =  +1  cyclotron  resonance  velocity  lies  in  regions 

dV 

^  3F 

of  large  anisotropy,  —  >  0. 

3a 

Figures  7  and  8  summarize  the  unstable  frequencies  and  wave  normal 
angles  for  the  electron  distribution  functions  shown  in  Figure  2.  Figure 
7  depicts  the  waves  that  are  generated  by  the  m  =  0  and  +1  resonances  for 
the  October  15,  1978,  velocity  distribution.  The  unstable  region  is 
indicated  by  the  shading.  These  waves  have  k||  <  0,  corresponding  to  wave 
vectors  directed  toward  the  magnetosheath.  The  unstable  frequencies  and 
wave  normal  angles  for  generation  via  the  m  =  -1  resonance  are  shown  in 
Figure  8.  In  the  bottom  panel,  the  results  for  the  October  15,  1978, 
distribution  are  shown  while  the  top  panels  correspond  to  the  two 
December  13,  1977,  distributions.  These  whistler  mode  waves  are  all 
generated  by  the  m  =  -1  cyclotron  resonance  term  and  have  k||  >  0. 

It  is  desirable  to  identify  the  growth  rates  that  are  large  enough 
to  produce  a  rapid  increase  in  the  intensity  of  whistler  mode  noise  at 
Earth's  bow  shock.  To  do  so,  the  critical  growth  rate,  wc,  required  to 
produce  large  growth  over  a  small  distance  is  calculated.  For  simpli¬ 
city,  the  critical  growth  rate  is  defined  to  be  the  growth  rate  necessary 
for  an  order  of  magnitude  increase  in  the  amplitudes  of  the  wave  electric 
and  magnetic  fields  over  a  distance  of  100  km.  An  order  of  magnitude 
increase  in  the  intensity  of  whistler  mode  noise  at  Earth's  bow  shock  is 
consistent  with  plasma  wave  observations  (e.g.,  see  Rodriguez  and  Gurnett 
[1975]).  The  spatial  length  over  which  the  growth  occurs,  100  km,  is  in 
agreement  with  bow  shock  scale  lengths  obtained  by  Russell  [1979]. 
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To  calculate  the  critical  growth  rate,  the  effect  of  the  solar 
wind  flow  on  the  propagation  of  a  whistler  mode  wave  packet  must  be 
taken  into  account.  This  is  because  the  solar  wind  speed  is  often 
comparable  to  the  whistler  mode  group  speed.  Because  only  a  rough 
estimate  of  the  critical  growth  rate  is  desired,  two  extreme  situations 
are  considered:  wave  packets  with  group  velocities  in  the  direction  of 
the  solar  wind  flow  and  wave  packets  with  group  velocities  in  the 
opposite  direction  of  the  solar  wind  flow.  The  critical  growth  rates 
for  these  two  cases  are  given  by 
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In  these  equations,  Vg  is  the  whistler  mode  group  speed,  vsw  is  the 
solar  wind  speed,  and  L  =  100  km  is  the  characteristic  scale  length  of 
the  transition  region.  The  group  speed  is  calculated  using 
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where  c  is  the  speed  of  light  and  0  is  the  wave  normal  angle.  The  solar 
wind  speed  is  taken  to  have  the  constant  value  400  km/s. 

Equations  5  and  6  are  used  to  calculate  the  critical  growth  rate  for 
waves  with  wave  vectors  directed  toward  the  magnetosheath,  k u  <  0,  and 
for  waves  with  wave  vectors  directed  toward  the  solar  wind,  k N  >  0.  The 
results  obtained  are  summarized  in  Figures  7  and  8.  The  region  contained 
inside  the  dashed  line  in  Figure  7  defines  the  frequencies  and  wave 
normal  angles  that  have  a  growth  rate,  a>j,  greater  than  the  critical 
growth  rate,  wc(k n  <  0).  For  wave  vectors  directed  toward  the  magneto¬ 
sheath,  k||  <  0,  waves  with  frequencies  between  about  5  and  85  degrees  and 
with  intermediate  wave  normal  angles  have  the  largest  growth  rates.  For 
wave  vectors  directed  toward  the  solar  wind,  k(|  >  0,  Figure  8  illustrates 
that  nearly  all  of  the  growth  rates  for  the  velocity  distributions 
measured  on  October  15,  1978,  and  at  1734:55.7  on  December  13,  1977,  are 
large.  However,  the  distribution  obtained  at  1734:58.5  on  December  13, 
1977,  does  not  contain  enough  free  energy  to  produce  growth  rates  larger 
than  the  critical  growth  rate.  At  this  time  it  is  not  known  whether  or 
not  this  is  a  common  feature  of  velocity  distributions  that  are  measured 
at  large  shock  penetration  lengths. 

It  should  be  noted  that  the  calculated  growth  rates  are  sometimes  so 
large  that  the  condition  uji/eDj.  <<  1  is  not  satisfied.  When  this  condi¬ 
tion  is  not  met,  the  growth  rate  calculated  using  Equation  1  is  not 
accurate.  It  is  expected  that  a  more  exact  calculation  of  the  growth 
rate  would  yield  a  smaller  positive  value.  While  this  difficulty  does 
point  out  the  fact  that  some  of  the  growth  rates  presented  in  this  paper 
are  rough  estimates,  the  signs  of  the  calculated  growth  rates  are 


accurate . 


2 


V.  DISCUSSION 

This  study  has  shown  that  electron  distribution  functions  measured 
using  the  fast  plasma  experiment  on  ISEE  2  within  the  transition  region 
of  Earth's  bow  shock  are  unstable  to  a  large  spectrum  of  whistler  mode 
waves.  The  unstable  whistler  mode  waves  have  plasma  rest  frame  frequen¬ 
cies  between  about  .1  and  100  Hz  and  have  wave  vector  directions  with 
respect  to  the  magnetic  field  which  range  from  0°  to  the  resonance  cone 
angle.  In  general,  the  lower  frequency  (.1  to  10  Hz)  waves  have  wave¬ 
lengths  between  about  500  and  100  km,  while  the  higher  frequency  waves 
(10  to  100  Hz)  have  wavelengths  between  about  100  and  20  km.  The  whistler 
mode  waves  are  generated  by  either  electromagnetic  Landau  resonance  or 

electromagnetic  cyclotron  resonance  or  both.  The  free  energy  regions  in 

k  ii 

the  electron  distribution  functions  are  regions  of  positive  -jy — r  — —  and 

1  k  II  '  '3V  I! 

regions  of  T^  >  Tfl  anisotropy,  both  attributable  to  the  field-aligned 
electron  beam.  If  it  is  assumed  that  a  background  continuum  of  whistler 
mode  noise  is  available  in  the  solar  wind  for  amplification,  as  indicated 
by  the  measurements  of  Neubauer  et  al.  [  1 9 7 7 J ,  the  stability  calculations 
indicate  that  many  of  the  unstable  frequencies  and  wave  normal  angles  have 
growth  rates  large  enough  to  give  rise  to  large  whistler  mode  intensities 
at  Earth’s  bow  shock. 

To  fully  evaluate  the  growth  of  whistler  mode  noise  in  the  how 


shock,  it  is  necessary  to  compute  the  growth  integrated  along  typical 
ray  paths  through  the  shock.  However,  a  detailed  treatment  of  this 


problem  is  difficult  because  the  magnetic  field  geometry  within  the 
transition  region  of  Earth's  bow  shock  is  complicated.  To  obtain  a 
qualitative  understanding  of  the  subsequent  propagation  of  the  generated 
whistler  mode  waves,  a  simple  model  for  the  shock  geometry  is  adopted. 
The  model  is  that  of  a  parallel  shock  with  the  upstream  solar  wind  velo¬ 
city  in  the  opposite  direction  of  the  shock  normal.  The  ramp  in  the 
magnetic  field  amplitude  defines  the  shock  transition  region  where  the 
electron  distribution  functions  shown  in  Figure  2  are  present.  The 
electron  distribution  function  measured  on  October  15,  1978,  represents 
the  distribution  function  often  observed  near  the  upstream  boundary  of 
the  magnetic  ramp.  The  electron  distribution  functions  measured  on 
December  13,  1977,  represent  the  distribution  functions  often  observed 
within  the  magnetic  ramp. 

The  waves  depicted  in  Figure  7  are  produced  near  the  upstream 
boundary  of  the  magnetic  ramp  by  the  m  =  0  and  the  m  =  +1  resonances  and 
have  k||  in  the  downstream  direction.  Because  the  waves  with  frequencies 
between  about  5  and  85  Hz  have  the  largest  growth  rates,  at  the  upstream 
shock  boundary  a  rapid  increase  in  whistler  mode  intensity  over  this 
frequency  range  is  expected.  For  the  parallel  shock  model  the  wave 
packets  propagate  downstream  and  are  convected  downstream  by  the  solar 
wind  flow.  This  is  due  to  the  fact  that  the  group  velocity  of  the 
whistler  mode  is  confined  to  angles  less  than  about  20°  from  the  mag¬ 
netic  field  [Helllwell,  1965],  As  these  wave  packets  move  downstream, 
they  will  encounter  distribution  functions  similar  to  those  measured  on 
December  13,  1977.  Because  the  electron  distribution  functions  measured 

on  December  13,  1977,  are  stable  to  waves  with  k ,,  <  ( > ,  the  waves  will  he 
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damped  as  they  move  downstream.  It  is  expected  that  most  of  these  waves 
will  not  reach  the  magnetosheath  as  the  damping  rate  can  be  large. 

The  waves  depicted  in  Figure  8  are  generated  by  the  m  =  -1  cyclo¬ 
tron  resonance  and  have  Vc n  in  the  upstream  direction.  For  the  simple 
shock  model  adopted  the  solar  wind  will  convect  the  wave  packets  down¬ 
stream  if  the  group  speed  is  smaller  than  the  solar  wind  speed,  which 
has  the  typical  value  400  km/s.  For  plasma  parameters  typical  of  the 
bow  shock  region  and  for  the  frequencies  of  interest  in  this  study,  the 
group  speed  of  the  whistler  mode  is  smaller  than  400  km/sec  for  wave 
normal  angles  near  0Res.  Consequently,  in  the  simple  shock  model  only 
the  whistler  mode  waves  with  wave  normal  angles  near  the  resonance  cone 
angle  are  convected  downstream  by  the  solar  wind  flow.  The  remaining 
whistler  mode  waves  either  travel  upstream  against  the  solar  wind  flow 
or  stand  in  the  shock.  It  is  possible  that  the  standing  whistlers  so 
produced  contribute  to  the  overall  shock  structure  because  many  of  the 
standing  whistlers  have  wavelengths  equal  to  the  shock  thickness. 

Because  the  electron  distribution  function  measured  on  October  15, 
1978,  represents  the  typical  distribution  near  the  upstream  shock 
boundary  the  growth  rate  calculations  for  this  distribution  predict 
which  whistlers  will  travel  into  the  solar  wind.  The  calculations 
summarized  in  the  bottom  panel  of  Figure  8  suggest  that  low  frequency 
(»7-30  Hz)  waves  travel  into  the  upstream  solar  wind.  A  range  in 
wave  normal  angles  from  0°  to  near  0ges  is  expected. 

These  predictions  are  summarized  as  follows: 

1)  Intense  broadband  whistler  mode  noise  is  generated  by  the 
m  =0,  +1,  and  -1  resonances  in  the  shock  transition 


region . 


2 


2)  The  whistler  mode  noise  in  front  of  the  magnetic  ramp  will 
be  most  prevalent  at  low  (7  -  30  Hz)  frequencies.  These 
waves  are  generated  by  the  normal,  m  *=  -1,  cyclotron  reso¬ 
nance,  and  have  wave  normal  angles  from  0°  to  near  0Res. 

3)  At  the  upstream  boundary  of  the  magnetic  ramp,  whistler  mode 
waves  with  frequencies  between  about  .1  and  100  Hz  are  pro¬ 
duced  via  the  m  =  0  and  +1  resonances.  The  waves  with  fre¬ 
quencies  between  about  5  and  85  Hz  have  the  largest  growth 
rates.  Because  these  waves  are  damped  as  they  travel  down¬ 
stream,  most  of  the  waves  are  not  expected  to  reach  the 
magnetosheath. 


It  is  desirable  to  compare  the  predictions  of  this  study  with 
observations  of  whistler  mode  waves  near  Earth's  bow  shock.  However, 
because  the  predicted  frequencies  are  with  respect  to  the  plasma  rest 


frame,  they  are  not  the  frequencies  actually  measured  by  a  plasma  wave 
receiver.  The  measured  frequencies  are  Doppler-shifted  by  the  relative 
motion  of  the  satellite  and  the  plasma.  Simple  calculations  show  that 
in  most  cases  the  measured  frequency  is  expected  to  be  between  about 
one-half  to  two  times  the  rest  frame  frequency.  Because  the  Doppler 
effect  at  the  shock  is  difficult  to  treat  in  detail,  in  this  paper  only  a 
qualitative  comparison  of  the  predictions  of  the  study  and  observations 
of  whistler  mode  waves  will  be  made. 

In  Figure  9,  data  obtained  on  December  13,  1977,  by  the  plasma  wave 
experiment  and  the  magnetic  field  experiment  on  board  ISEE  1  are  shown. 
ISEE  1,  as  opposed  to  ISEE  2,  data  are  shown  because  wave  magnetic  field 
measurements,  which  are  most  useful  for  detecting  whistler  mode  noise, 


were  only  available  from  ISEE  1  for  the  event  studied.  The  spatial 
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separation  of  the  two  satellites  is  small.  By  comparing  the  ISEE  1  and 
the  ISEE  2  magnetic  field  profiles  at  these  times,  the  temporal  separa¬ 
tion  of  the  two  data  sets  measured  on  December  13,  1977,  is  found  to  be 
about  15  seconds.  For  example,  the  time  1735:00.0  on  ISEE  2  corre¬ 
sponds  to  1735:15.0  on  ISEE  1. 

The  top  panel  of  Figure  9  shows  the  wave  magnetic  fields  in  ten 
channels  covering  the  frequency  range  5.6  Hz  to  1  kHz  and  the  time  period 
1730:00.0  to  1740:00.0  UT  on  December  13,  1977.  The  bottom  panel  depicts 
the  magnetic  field  strength  over  the  same  time  period.  The  magnetic  ramp 
defining  the  shock  transition  region  is  clearly  visible,  with  the 
upstream  boundary  at  about  1735:05.0.  At  this  time,  the  satellite  enters 
the  shock  from  the  solar  wind.  It  is  assumed  that  the  electron  distribu¬ 
tion  function  at  the  upstream  boundary  of  the  magnetic  ramp  is  similar  to 
the  October  15,  1978,  electron  distribution  function.  Within  the  mag¬ 
netic  ramp  at  about  1735:20.0,  the  electron  distribution  functions  are 
those  measured  by  the  ISEE  2  fast  plasma  experiment  on  December  13,  1977. 

In  Figure  9,  whistler  mode  noise  is  visible  in  the  bottom  six  fre¬ 
quency  channels  of  the  wave  magnetic  field  data.  As  predicted  by  this 
study,  the  noise  has  largest  intensity  within  the  magnetic  ramp.  The 
frequency  range  over  which  the  noise  is  observed,  from  about  5.6  to  100 
Hz,  is  in  good  agreement  with  the  results  of  this  study.  In  the  solar 
wind  immediately  upstream  of  the  magnetic  ramp,  the  noise  is  most 
intense  at  low  frequencies,  in  agreement  with  prediction  2  above. 

This  figure,  which  is  typical  of  whistler  mode  observations  at  the  bow 
shock,  illustrates  that  the  predictions  of  this  study  and  the  observa¬ 
tions  of  whistler  mode  waves  compare  favorably. 
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Figure  1 


Figure  2 


Figure  3 


FIGURE  CAPTIONS 


This  figure,  reproduced  from  Feldman  et  al.  [1982], 
shows  parallel  slices  of  electron  distribution  func¬ 
tion  data  measured  by  the  ISEE  2  fast  plasma  experi¬ 
ment.  As  the  satellite  moves  from  the  solar  wind 
into  the  magnetosheath,  field-aligned  electron  beams 
are  observed. 

This  figure  depicts  contour  plots  of  the  logarithm  of 
the  analytic  fits  to  the  electron  distribution  func¬ 
tion  data  analyzed  in  this  study.  These  distribu¬ 
tion  functions  were  measured  within  the  transition 
region  of  Earth's  bow  shock.  A  modified  Lorentzian 
models  the  flat-topped  component  of  the  electron  dis¬ 
tribution  and  a  convected  Maxwellian  with  a  T^  >  T n 
anisotropy  models  the  field-aligned  electron  beam. 

The  left  column  of  this  figure  presents  a  sample  of 
the  growth  rates  obtained  for  the  electron 
distribution  function  measured  on  December  13,  1977, 
at  1735:55.7.  Plotted  is  the  ratio  of  the  total 
growth  rate  to  real  frequency.  These  waves  are 
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Figure  4 


Figure  5 


generated  by  the  m  =  -1  cyclotron  resonance  and  have 
k  ||  >  0.  The  right  column  of  this  figure  depicts  the 
m  =  -1  cyclotron  resonance  velocity  as  a  function  of 
wave  normal  angle  for  the  two  frequencies  of 
interest . 

In  this  figure,  the  range  in  the  m  =  -1  cyclotron 
resonance  velocity  corresponding  to  instability  for 
u>t/2tt  =  110  Hz,  as  obtained  from  Figure  2,  is  over- 
layed  on  a  contour  plot  of  the  electron  distribution 
function  measured  on  December  13,  1977,  at  1734:55.7. 
The  m  =  -1  cyclotron  resonance  velocities  lie  in 
regions  of  large  anisotropy,  — —  <  0. 

oCC 

The  left  column  of  this  figure  presents  a  sample  of 
the  growth  rates  obtained  for  the  electron  distribu¬ 
tion  function  measured  on  October  15,  1978,  and  for 
the  frequency  ^r / 2 tt  =  80  Hz.  The  unstable  waves  have 
k  u  <  0,  corresponding  to  wave  vectors  directed  toward 
the  magnetosheath.  The  top  left  panel  shows  the  m  =  0 
Landau  contribution  to  the  growth  rate  while  the 
bottom  left  panel  shows  the  m  =  +1  cyclotron  contri¬ 
bution  to  the  growth  rate.  For  this  frequency  and 
these  wave  normal  angles  the  m  =  -1  cyclotron  contri¬ 
bution  to  ^i/^r  are  negative  and  negligible.  The 
right  column  depicts  both  the  m  =  0  Landau  resonance 
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Figure  6 


Figure  7 


Figure  8 


velocity  and  the  m  =  +1  cyclotron  resonance  velocity 
as  a  function  of  the  wave  normal  angle. 

In  this  figure,  the  range  in  the  m  =  0  Landau  and  the 

m  =  +1  cyclotron  resonance  velocities  corresponding 

to  instability  for  the  case  ^4* / 2 tt  =  80  Hz,  as 

obtained  from  Figure  4,  is  overlayed  on  a  contour 

plot  of  the  electron  distribution  function  measured 

on  October  15,  1978.  The  m  =  0  Landau  resonance 

k  !l  3F 

velocities  lie  in  regions  of  positive  ■. -  1  — • —  while 

|  k  ||  |  3v  n 

the  m  =  +1  cyclotron  resonance  velocities  lie  in  a 

3F 

region  of  large  anisotropy,  —  >  0. 

da 

This  figure  depicts  the  unstable  frequencies  and  wave 
normal  angles  that  are  generated  by  the  m  *  0  and  +1 
resonances  for  the  electron  distribution  function 
measured  on  October  15,  1978.  The  waves  have  k n  <  0, 
corresponding  to  wave  vectors  directed  toward  the 
magnetosheath.  The  largest.  ferowth  rates  lie  in  the 
shaded  region  bounded  by  the  dashed  line. 

This  figure  depicts  the  unstable  frequencies  and  wave 
normal  angles  that  are  generated  by  the  m  =  -1  cyclo¬ 
tron  resonance  for  the  electron  distribution  functions 
shown  in  Figure  2.  All  of  the  waves  and  have  k ,,  >  0, 
corresponding  to  wave  vectors  directed  toward  the 


solar  wind.  The  growth  rates  for  the  distribution 
measured  at  17  34:58.5  on  December  13,  19/7,  art- 
smaller  than  the  critical  growth  rate. 

This  figure  depicts  data  obtained  with  the  ISKK-I 
plasma  wave  experiment  and  the  magnetic  field  experi¬ 
ment.  The  top  panel  shows  the  wave  magnetic  fields 
for  six  frequency  channels.  The  bottom  panel 
illustrates  the  magnetometer  data.  The  magnetic  ran} 
at  the  shock  is  clearly  visible.  The  characteristics 
of  the  whistler  mode  nc.se  at  the  shock  are  in  good 
agreement  with  the  predictions  of  this  study. 
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